Majorana Fermions and Odd-frequency Cooper Pairs in a Nano Wire 
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We discuss a strong relationship between Majorana fermions and odd-frequency Cooper pairs 
which appear at a disordered normal (N) nano wire attached to a topologically nontrivial supercon- 
ducting (S) one. The transport properties in superconducting nano wire junctions show universal 
behaviors irrespective of the degree of disorder: the quantized zero-bias differential conductance at 
2e^/h in NS junction and the fractional current-phase {J — ip) relationship of the Josephson effect in 
SNS junction J oc sin(i/j/2). Such behaviors are exactly the same as those found in the anomalous 
proximity effect of odd-parity spin-triplet superconductors. We show that the odd-frequency pairs 
exist wherever the Majorana fermions stay. 

PACS numbers: 74.25.F-, 74.45. +c, 74.78. Na 
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Majorana fermion (MF) satisfying a special relation 
of 7 = 7twith 7 (7^) being the annihilation (creation) 
operator has been an exciting object since the original 
prediction by Majorana-. Finding of MFs and control- 
hng of Majorana bound states (MBSs) are hot research 
issues in condensed matter physics^ from the view of 
potential application of MBS to the topological quan- 
tum computationally. To date, we have known several 
promising systems hosting MFs such as spin-triplet p- 
wave superconductors^"*, topological insulator /super- 
conductor heterostructures^emiconductor / supercon- 
ductor junctions with strong spin-orbit coupling iSr^, 
helical superconductorsi^, and superconducting topolog- 
ical insulators-?. Most attracting case among them is 
the semiconductor nano wire fabricated on top of a su- 
perconductor because of its easy tunability of MBS by 
changing the chemical potential in the nano wire and by 
applying the Zeeman field onto it^i^iil. Actually, a plenty 
of theoretical studies have discussed MFs or MBSs in 
such nano wire a-'^^d^ . The zero-bias conductance peak 
reported very recently would be considered as an ev- 
idence of MFs (MBS) i&iiliSa. These researches have 
stimulated a number of theoretical investigation on un- 
usual charge transport phenomena through the MBS in 
normal-metal / superconductor (NS) and superconduc- 
tor /normal- metal/ superconductor (SNS) junctions on 
nano wires^. However, no studies have ever tried to an- 
alyze features of Cooper pairs which support the anoma- 
lous transport properties. We address this issue in the 
present Letter. 

Odd-frequency Cooper pairing was originally proposed 
to understand nature of unconventional superfluidity and 
superconductivitjj2^. Ubiquitous appearance of the odd- 
frequency pairs at the surface of superconductors and 
near the interface of superconducting junctions has been 
established and widely accepted in recent years^^i^i. The 
zero-energy Andreev bound state (ABS) at the surface 
of unconventional superconductors ^® is reinterpreted in 
terms of the odd-frequency Cooper pairing^^. In particu- 
lar, the odd-frequency Cooper pairs make the background 



of the anomalous proximity effect in a diffusive normal 
metal attached to a spin-triplet superconductor^^: (i) 
the large zero-energy quasiparticle density of states in 
a normal metaP®''^^, (ii) the quantized zero-bias conduc- 
tance at twice of the Sharvin's value in diffusive NS junc- 
tions—, (iii) the fractional current( J)-phase((/3) relation- 
ship of J oc sin((p/2) in diffusive SNS junction^, (iv) the 
zero-bias anomaly in nonlocal conductance spectr and 
(v) the anomalous surface impedance in NS bilayers'^^. 

In this Letter, we show that disordered NS and SNS 
junctions of nano wire indicate the properties of (i)-(iii) 
when the superconducting nano wire is topologically non- 
trivial. In addition, the amplitude of odd-frequency pairs 
in the normal nano wire suddenly grows as soon as the 
superconducting nano wire undergoes the transition to 
topologically nontrivial phase. The unusual transport 
phenomena due to the MFs2S are nothing but the anoma- 
lous proximity effect due to the odd-frequency pairs. We 
will conclude that the odd-frequency Cooper pairs are 
indispensable to realizing MFs in solids. 
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FIG. 1: (color online), 
junctions. 



(b) SNS nano wire 



Schematic pictures of NS and SNS 



Let us consider a nano wire with strong spin-orbit cou- 
pling fabricated on a junction of an insulator and a metal- 
lic superconductor as shown in Fig. [TJ A segment on the 
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insulator and that on the superconductor are in the nor- 
mal and the superconducting states, respectively. The 
diameter of the nano wire is sufficiently small so that 
the number of propagating channel is unity for each spin 
degree of freedom. We describe the present nano wire 
by using the tight-binding model in one-dimension, for 
noninteracting electrons^'^, 
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l<j<L,a 

i/..=^[Ae»-^44, + H.c. 
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where c] ai'^j^a) is the creation (annihilation) operator of 
an electron at the lattice site j with spin a = (f or J,), t 
denotes the hopping integral, /i is the chemical potential, 
and A is the pair potential in the superconducting seg- 
ment. The Pauli matrices in spin space are denoted by aj 
ioT j — 1 — 3 and the unit matrix of 2 x 2 is ctq. The on-site 
potential in the normal segment is given randomly in the 
range of —W/2 < Vj < W/2. We measure the energy and 
the length in units of t and the lattice constant, respec- 
tively. Throughout this paper, we fix several parameters 
a.s fi = t, W — 2t, and the pair potential at the zero 
temperature A = O.Olt. The number of samples used 
for the random ensemble averaging is typically 10'^ — 10^. 
By tuning the magnetic field B as shown in Fig. [Tl it is 
possible to introduce external Zeeman potential Vex- For 
Vex > Vc= \/A§+/?, the number of propagating chan- 
nels becomes unity and the superconducting segment un- 
dergoes the transition to topologically nontrivial phase. 
In the tight-binding model, the finite band width gives 
an additional condition Vex < Vc2 = 4t — /i for the topo- 
logical phase. Here we briefly summarize calculated re- 
sults of the normal conductance of the disordered nano 
wires with using the recursive Green function method^. 
By analyzing the Hamiltonian Hq + Hd, we confirmed 
that the normal conductance decays exponentially with 
increasing iM-. This is because one-dimensional disor- 
dered wires are basically in the localization regime. 

At first, we calculate the differential conductance Gns 
of NS junctions based on the standard formula—, 
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where we consider the Hamiltonian Hq + Hd + Hg. In 




Eq. ©, 



and 
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are the normal and Andreev re- 
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FIG. 2: (color online). The differential conductance of NS 
nano wires is plotted as a function of the bias-voltage for 
non-topological nano wire {Vex = A = 0) in (a) and for the 
topological nano wire ( Vex = 1.5f and A = 0.5t) in (b). In 
(b), we also plot the results for W/t — 4 and L = 10 with 
a solid line. The local density of states at the center of the 
disordered segment {j = 10) are shown for several Vex in (c) 
with A = 0.5t, where we add a small imaginary part iSe with 
Se — O.OIA to the energy. 



flection coefficients of the junction at energy E mea- 
sured from the fermi level. We show Gns in units of 
Gq — 2e^/h as a function of the bias- voltage eV for 
the non-topological and the topological nano wires in 
Figs, mja) and (b), respectively. The length of disor- 
dered segment L is chosen as 10, 20, and 50 lattice con- 
stants. The conductance for the non-topological nano 
wires ( Vex = A = ) in (a) decreases with increasing 
the length of disordered segment L for all eV. The sim- 
ilar tendency can be seen also in the conductance of the 
topological nano wires ( Vex = l-5i and A = 0.5t) in (b) 
for finite eV. However the zero-bias conductance of the 
topological nano wires is quantized at Gq irrespective of 
L, which is an intrinsic phenomenon in the presence of 
MF.The results suggest a perfect transmission channel 
due to the penetration of a resonant state into the disor- 
dered segment. The local density of states (LDOS) in the 
disordered nano wire supports this statement as shown 
in Fig. [21(c), where we plot the LDOS at the center of 
the disordered segment (j = 10) as a function of E for 
L = 20 and A = 0.5t. The results are normalized to the 
density of states at the fermi level in clean normal nano 
wire iVo. The LDOS for Vex > Vc show the large zero- 
energy peak reflecting the MBS as shown in the results 
for Vex/t = 1.5 and 1.2. On the other hand, the LDOS 
for Vex/t ~ 0.5 and 0.8 are almost flat around the zero 
energy. 

Secondly, we explain why the superconducting nano 
wire shows the anomalous proximity effect which is 
unique to the px-wstve spin-triplet superconductor. The 
single particle Hamiltonian in Eq. ([l} is essentially equiv- 
alent to ft-o(^) = 'CfcO'o — Xka2 — VexO'3 in momentum space 
with — h^k^ /{2m) — /i. By applying a unitary trasfor- 
mation diagonalizing /lo(^) smd — /i5(— fc), Bogoliubov-de 
Gennes (BdG) Hamiltonian of the nanowire reduces to 
2x2 Hamiltonian for Vex > Vc (See also Appendix A for 
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details), 



ho{k) «Ao-2 



6 - A Afee"/2 



, (6) 



with A = yV^IT+T^^O^ A^ = AXk/A because a spin 
branch pinches off from the fermi level (i.e., ^ + ^ > 0). 
The right-hand side of Eq. © is equivalent to the BdG 
Hamiltonian of spin less p^j-wave superconductor in one- 
dimension when we focus on low energy excitation after 
redefining the chemical potential fi + A ^ fi and the 
pair potential Afce*''/^ A{k/kF) with kp being the 
fcrmi wave number. Therefore physics in the topological 
nanowire is the same as that of p^j-wave superconductor. 
In fact, we have confirmed that the Josephson current in 
SNS junctions of disordered nanowire show the fractional 
current-phase relationship at low temperature as shown 
in Fig. |3lja)--. Here we attach the second supercondutor 
for j < by adding Hs2 to Hq + Hg + Hd, and plot the 
Josephson current J as a function of Atp — Lp2 ~ aX 
T = O.OOlTc for Vex = IM and A = 0.5t. The results 
show J oc sin(A(/?/2) for — tt < Atp < tt irrespective of 
L. For comparison, we also plot the results for Vex = 
A = and L = 50 with a solid line which shows usual 
sinusoidal current phase relationship. Correspondingly 
the Josephson critical current plotted as a function of 
temperature shows so called low temperature anomaly22- 
in Fig.[3i;b). 

Next we discuss the relationship between Majorana 
fermions and odd-frequency Cooper pairs by analyzing 
the Green functions in junctions of p^j-wave supercon- 
dutor. A semi-infinite wire of p^^-wave superconduc- 
tor occupying x > Q hosts a Majorana fermion around 
its edge a.i E — 0. Solving the BdG equation, the 
wave function of such surface state is calculated to be 
4>o[xY = [u^[x),vo {x)]^ , where uo{x) = C{x)x, vo{x) = 
Cix)x*, C(x) = y27^e-^/2««sin(fcx), x = e'^'/'^e'^^^ , 
and ^0 is the coherence length. The electron operator 
includes the contribution from such surface state ipo{x) 
as represented by -00 (2;) = Xli^): V'oC^^) — X*'l'{x) with 
7(0;) = C{x){'~fo + 7o)- Here 70 is the annihilation oper- 
ator of the Majorana bound state. The special relation 
Vo{x) — Uq{x) plays a crucial role in the Majorana re- 
lation of 7(x) = 7^(0;)^'^. As a result, the two Green 
functions calculated for \E\ ^ A 

g[-Cj, X,X ) ~ zjii-r ' V/ 

L + lOe 

f(j^ ,^ _ Uf){x)vl{x') + vI{x)uq{x') 
J(Jij,X,X , [6 

h + lOc 

depend on each other. Since vq{x) — Uq{x), they satisfy 

g{E, X, x') = {x*ff[E, X, x') = I{E- x, x'). (9) 

This relation directly bridges Majorana fermions and 
odd-frequency Cooper pairs. The real (imaginary) part 
of {x*)'^f{E;x,x) is an odd (even) function of E, which 



represents the odd-frequency symmetry of Cooper pairs. 
The orbital part is s-wave symmetry when / is calcu- 
lated at a; = x' . In fact, the imaginary part of g{E; x, x) 
must be even function of E because it represents LDOS 
of the Majorana bound state. It is possible to check this 
argument in a junction which consists of a normal metal 
(x < 0) and a pa;-wave superconductor {x > 0) in one- 
dimension'^^- (See also Appendix B). At the NS interface 
{x = 0), we introduce a potential barrier VoS{x) whose 
normal transmission coefficient is t„ — kp/{kp + izq) 
with zq — Vq/Tivf- When we focus on the subgap energy 
\E\ <C A in the tunneling limit |t„| ^ 1, we find that 
the Green functions in the superconductor x > Q satisfy 
Eq. ([9]) and become 



I{E] X, x) 



ttNqA 



E + iA\tn\y2 



e-^/«0sin2(fcx). (10) 



For comparison, we show the anomalous Green function 
in a uniform ^^.-wave superconductor 



{x*ff{E;x,x')=-t 



, ttNo a sin k{x — x') 



VWT70+)2~A2' 



(11) 



with ^[E + iQ+Y - A2 being sgn{E)^/E'^ - A^ for 
\E\ > A and iVA^ — E^ for \E\ < A, where we assume 
|a; — x'l ^ ^0- The anomalous Green function satisfies 
f{x — x') = —f{x' — x) reflecting the odd-parity symme- 
try. In contrast to Eq. (fTO)) . the real (imaginary) part 
of ix*)"^ f{E, X, x') is an even (odd) function of E, which 
represents the even-frequency symmetry. 

The important relation in Eq. Q can be confirmed in 
the normal segment of NS nano wire as shown in Fig. |4l 
where we fix the energy at -E = and plot gx\ij,j) and 
~/ttO' j) tli^ center of the normal segment j = 10. We 
note that an extra phase factor (p — 7t/2 in Eq. © makes 
(X*)'^ = -1 in Eq. For Vex > the results show 
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FIG. 3: (color online), (a): Current-phase relationship in SNS 
junctions of topological wire at T/Tc = 0.001 for Vex = 1.5t 
and A = 0.5t. For comparison, the results for non topological 
wire (Vex = A = and L — 50) is plotted with a solid line, 
(b): The Josephson critical current versus temperature in a 
topological nano wire. 

Im((7||) = — Im(/||-). The real part of /^j is always zero at 
E = due to the odd-frequency symmetry. Correspond- 
ingly, Re(gii) also goes to zero for Vex > Ve- In addition. 
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— Im((/i-f) = Im(/||) suddenly increases as Vex increasing 
across Vc, which corresponds to the zero-energy peak in 
LDOS in Fig. EJc). The results demonstrate the pene- 
tration of Majorana fermions and odd-frequency Cooper 
pairs into the normal disordered segment at the same 
time. 
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FIG. 4: (color online). The normal gy^ and the anomalous 
/ti Green function for _E = at j = 10, with L = 20. We 
introduce a small imaginary part 5^ = O.OOIA. 

In the NS junction of p^j-wave superconductor, it is 
possible to calculate exactly the wave function in the 
presence of a single impurity Vi5{x — Xi) in the normal 
metal by solving the Lippmann-Schwinger equation^ 
(See also Appendix C), 



(t)n{.x) =(t)^;:^Hx) + G{E;x,x{)Viaz(l)n{x{), 



(12) 



where (/)„ {x) is the wave function in the presence of the 
impurity and (x)]"^ 



that in the ballistic case with r^e and The being the nor- 
mal and the Andreev reflection coefficients at the NS in- 
terface from the electron branch, respectively. By putting 
X = Xi, the equation in Eq. (|12p has the closed form for 
4>n{xi), which results in 



0„(a:) = C^(x) 



1 

Y 



Zi{zi + i)X 



e**"(l-r)r^e 



, (13) 



at £' = for x < Xi, where Y = 1 + zfX, X = I — rhereh, 
Zi = Vi/hvp, and r^h is the Andreev reflection coefficient 
of the NS interface from the hole branch. We have al- 
ready taken into account the absence of the normal reflec- 
tion at the NS interface at E — (i.e., Tee = rhh = 0). At 
E — 0, the Andreev reflection coefficients r^e = — ie^*"^ 



and Tp 



satisfy an important relation VehThe = 1 



Teee *'=^,r;iee'''^]'^ IS 



which eliminates the second term of Eq. (I13|) . Thus the 
zero-bias conductance quantization at Gq = jh holds 
even in the presence of an impurity. The relation of 
TehThe = 1 IS nothing but the condition for forming the 
MBS at = 0. To discuss whole effects of scatterings by 
many impurities, we need to solve the nonlinear quasi- 
classical Usadel equation^Ss^E. The analytical expression 
of the zero-bias conductance also show Gi^s — Gq (See 
Appendx C). The diffusive normal metal is assumed in 
the Usadel equation. The validity of Gns = Gq in the 
localization regime is confirmed in numerical calculation 
in Fig. m^a) and in Fig. |S] 

In summary, we have theoretically discussed the 
anomalous transport phenomena in NS and SNS junc- 
tions of nano wires in which the superconducting seg- 
ment is topologically nontrivial and the normal segment 
is disordered by random impurity potential. The physics 
behind the anomalous transport can be understood in 
terms of the odd-frequency Cooper pairing. We conclude 
that Majorana fermions and odd-frequency Cooper pairs 
in solids are two sides of a same coin. 

This work was supported by KAKENHI on Innovative 
Areas "Topological Quantum Phenomena" and KAK- 
ENHI(22540355,22103005) from MEXT of Japan. 



Appendix A: Transformation of Hamiltonian 



The starting Hamiltonian of this paper is equivalent to 



NW 



2 1,2 



2m 



(Al) 
(A2) 



where is the chemical potential. Vex is the Zeeman potential due to external magnetic field, Xka2 represents the 
spin-orbit coupling, (Tq is the unit matrix in spin space, and aj for j = 1 — 3 are the Pauli's matrices. By applying 
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following unitary transformation'^^, the Hamiltonian is deformed as 

D^^HnwD = 







/^^e^{v+^/2) Aye*^ 

-Aye'^^ -^k + A 



D = 



Xk 

1' 


Av = 


A-^ 

A 


u ■ 
u 




U = 










)• 



a isgn{k)l3 
i sgn(/c)/3 a 



(A3) 

(A4) 
(A5) 

(A6) 



The topologically nontrivial phase is characterized by Vex > V/P + A^. In such case, only one dispersion remains at 
the fermi level for each Nambu space (i.e., ^k + A > Q). Therefore the Hamiltonian reduces to 2 x 2 Nambu space as 



H 



NW2 



^k-A AAe*(^+''/2) 



(A7) 



When we focus on the low-energy quasiparticle excitation, this is equivalent to the Hamiltonian describing the equal 
spin-triplet (spin less) p^^-wave superconductor 
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(A8) 



Here we redefine fj, + A ^ in the diagonal term and Axe^^l"^ Aikjkp) in the off-diagonal term. In previous 
papers2ii21i^i22ri2^, we have studied the anomalous proximity effect starting from the Hamiltonian in Eq. (|A8[) . 



Appendix B: Analysis of p^^-wave superconductor 
1. Green function and its representation 

The retarded Green functions are defined by the standard way 

G{x,t;x't') F{x,t;x't') 
Flx,t;x't') G{x,t-x't') 



(Bl) 
(B2) 



where ^/'(x)(V'^(x)) is the annihilation (creation) operator of a spin less electron. In the case of spin-triplet supercon- 
ductors, the electron operators are represented by the Bogoliubov transformation 



i/>(x) 
i^\x) 



E 



Ui,{x) vl(x) 



(B3) 



where is the annihilation operators of Bogoliubov quasiparticle belonging to Ei,. The wave function Uy{x) and 
v^{x) are obtained by solving the Bogoliubov-de Gennes (BdG) equation. The Green functions are expressed in 
spectral representation as 



G{E-x,x')=Y^ 
F{E;x,x')^Yl 



u^{x)ul{x') v*{x)v^{x') 



E + iS- E^ E + iS + E^ 
Uv{x)vl{x') v*{x)u^{x') 



E + iS- E^ E + iS + E^ 



(B4) 
(B5) 



where iS is a small imaginary part. 
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2. Uniform superconductor 



The retarded Green function of a uniform spin less p^^-wave superconductor is calculated to be 



E + ^ Asx 
Asx E -fl 



E-Q -Asx 
-Asx E + Q 



-ik \x—x'\ 



r(F.^ ^'^-( G{E;x,x') F{E-x,x')\ 
(.\,i^,x,x ) - ^ p^-^.^^^,^ G{E;x,x') J ' 

= k(l±^^ , n = ./{E + i5f^^, $ = diag[e'^/2^ 6-'"^'% = sgn{x - x'), 



(B6) 
(B7) 
(B8) 



where A^o is the density of states (DOS) at the fermi level in the normal state. The Green functions are calculated as 



G{E;x,x') 

-ie-''^F{E;x,x') 



.TtNq ikQ\x-x'\/{2ix) 



E 

— cos A;(x — x') + isin A;|a; — x'\ 



From the normal Green functions, the local density of states (LDOS) is calculated to be 

N{E,x) = —ImTrG{E,x,x) = NoRe^. 
The LDOS is an even function of E because of the relation 



_ ^2 A < 

^/{ETW^^ = { W A^ - E^ < < A . 

WE"^ - A2 E<-A 



(B9) 
(BIO) 

(Bll) 
(B12) 



It is evident that there is no subgap state in uniform superconductor. From the off-diagonal part, it is possible to 
check the pairing symmetry. The anomalous Green function satisfies F{x — x') = —F{x' — x), which indicates the 
odd-parity symmetry. In addition, the real part of —ie~'"^F{E, x, x') is an even function of E, whereas the imaginary 
part of it is an odd function of E. This means that Cooper pairs have the even-frequency symmetry. 



3. Majorana surface bound state 



Next we consider a semi- finite Pa;-wave superconductor which occupies x > as shown in Fig. 1(a). By solving the 
Bogoliubov-de Gennes equation, the wave function for subgap state is expressed by 



E + m 

Ae-^'P 



E-iVL 

-Ae-"^ 



g-ifcXg-x/2{o 



(B13) 



where = vA^-E^, A and B are constant. From the boundary condition at x = (i.e., '^s{x = 0) = 0), we find 
that a subgap state exists at i? = and that the wave function of it becomes 



uo{x) 
vo{x) 



= C{x) 



gi7r/4gi(^/2 
g-i7r/4g-i¥'/2 



C{x) = 



sin(fca;) 



We note that two components in the wave function satisfy an important relation 

uo{x) = vl{x). 

The BdG transformation reads, 



il){x) 

4>\x) 

(j>o(x) 
cl>lix) 



u^{x) v;{x) 

v^{x) ul{x) 



E 

Uo{x) v^{x) 
Vo{x) u^{x) 



^q{x) 



<l^l{x) 



70 
7^ 



(B14) 

(B15) 

(B16) 
(B17) 
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FIG. 5: System under consideration, (a): a semi-infinite p^-wave superconductor, (b): a clean normal- metal/superconductor 
(NS) junction of pa;-wave symmetry, (c): an impurity is introduced in the normal metal, (d): a diffusive normal metal is 
introduced in the NS junction. 



where 70 is the annihilation operator of the bound state. Together with Eqs. (IB14|) and (jB15|) . we find 

M^) =e*"/4e*^/2^(a;), (t>l{x) = e-'^'^^e-'^/^-fix), (B18) 
7(x) =C(x)(7o+7^). (B19) 

The fermion operator j{x) satisfies the Majorana relation 7(2;) = 7^(2;). When we focus on \E\ <C A, the contributions 
from V — {) become dominant in Eqs. (|B4p and (jB5|) . Near i? = 0, the normal and the anomalous Green functions 
satisfy a relation 

G{E; X, x') = -ie-'fF{E; x, x'), (B20) 
because they are calculated from Eq. ()B17P as 

L + 10 L + Id 

Eq. (IB20)) directly relates Maj orana fermions and odd-frequency Cooper pairs. When we consider x = x' , F{E,x,x) 
represents the pairing function of s-wave symmetry. The real part of —ie~^^F{E] x, x) is an odd-function of E and 
the imaginary part of it is an even function of E. This indicates that Cooper pairs have the odd-frequency symmetry. 
It is possible to confirm Eq. (|B20[) in NS junctions as discussed below. 



4. NS junction of superconductor 



It is possible to calculate the Green function of a junction which consists of a normal metal (x < 0) and a px-wave 
superconductor (a; > 0) in one dimension as shown in Fig. 1(b). In the superconductor, the retarded Green function 
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becomes^ 



GssiE;x,x') = $z 



U2 


UVSx 


UVSx 




1 —uv 






—uv 


} - 


~uv \ 


\ uv - 





NqttE 

ik'^ \x — x' 



V -UVSx \ ^_ik-\a:-x'\ 

— UVSx 



) —ik x+ik^ x' he , [ ^ | ik^ x—ik x') eh 

^ \ V^ UV J 



uv u 



^2 \ e~'''"(^+=^'V'''' 



(B22) 



'1 



+ r! = V(^^ + *'5)'-A2, fc± = fc(^l±^j, =sgn(x-x') (B24) 



17 \ 



(B23) 



for a;, x' > 0. The subscript ss in the Green function Gss{E\ x, x') means that both a; > and x' > indicate places 
in the superconductor. The normal and Andreev reflection coefficients are given by 



d^=^(2-it„n 



SS 1 ss 

k 



-Xu^-v^), r 



hh 



-{U^'V^ 



k + izo ' 



: — , Zo = Vo/fiVF 



k + izQ 



(B25) 
(B26) 



where tn and r„ are the normal transmission coefficients due to the potential barrier at the interface described by 
Vod{x). When we focus on the subgap energy in the tunneling limit (i.e., ^ A and |t„| ^ 1), we find, 



G,siE;x,x) = -ie~''^F,siE;x,x) « ttTVq 



A 



E + iA\tn\y2 



(B27) 



The imaginary part of Gss{E] x, x) gives the local density of the Andreev bound state and must be an even function 
of E. Therefore the real part of —ie~^^Fss{E]x,x) is an odd-function of E and the imaginary part of it is an even 
function of E, which indicates the odd-frequency symmetry. The condition in Eq. (jBlSp leads to the Majorana relation 
in operators at Eq. (jBlSp . Eq. (jB18l) results in Eq. (jB20l) . Then Eq. (|B20[) guarantees the odd- frequency symmetry 
of Cooper pairs. The orbital part of Cooper pairs is s-wave symmetry because —ie^'^^'^Fss is calculated aX x — x' . 



Appendix C: Perfect transmission at i5 = 

In the normal metal of NS junction < 0), the Green function is given by^ 



Gnn (^E ^ X, X ) 



-ik\x—x'\ I ^—ik(x-\-x' )^ee 



^ — ikx ~\~ikx /yt^fx 
^ ' nn 

^ik\x—x'\ _|_ ^ik{x-\-x') ^hh 



hh n 

nn r— . ; 



Gnji , X, X ) Ffiji , X, X ) 
Fjini^-^ 1 X, X ) Gjijii^E^ X, X ) 



(CI) 

(C2) 

(C3) 
(C4) 



The subscript nn in the Green function G„„(£'; x, x') means that both x < and x' < indicate places in the normal 
metal. We also confirmed the relation between the Green function in the normal metal, 



Gnn{E; X, x) = -ie "■'^F^niE; x, x) « TriVo 



A|t„| 



i? + iA|i„|2/2' 



(C5) 
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for \E\ <C A and |t„p <C 1. As we discussed in Eq. (IB20I) . this suggests the presence of Majorana fermions and 
odd-frequency Cooper pairs in the normal metal. 

It is possible to calculate exactly the wave function where a single impurity Vi6{x — Xi) exist in the normal metal 
as shown in Fig. 1(c) by using the Lippmann-Schwinger equation, 



(j)n{x) =(?!)™'(a;) + Gn„{E; X, Xi)Via3(l)n{xi), 

lini / \ 



ikx he 



(C6) 
(C7) 



where 0^"*(a;) is the wave function in the ballistic case and (j>n{x) is that in the presence of the impurity. By solving 
this equation at x = x^, we obtain 



The wave function for x < Xi in the presence of the single impurity is expressed by 

r . -1 -1 . . 

<j)n{x) =(f>'^'{x)+Gnn{E;X,X,)Via3 1- Gnn{E■,X^,X,)V^a3 (t^Ti^i)^ 



Y 



Y =l + z,{e 



i2kxi „hh 
nn 



i2kxi ee \ i ^2 



he 
nn 



he eh 
nn nn 



he eh , i2kxi hh . ~i2kxi ee \ 
nn nn ' nn ~^ nn ) ' 



nn/ ' V-^ 



Ai E — Q, the reflection coefficients become 



y,he 



-le 



-lip 



r''^ = 0. 

nn ^ 



These relations immediately lead to 



(C8) 

(C9) 

(CIO) 

(Cll) 
(C12) 

(C13) 

(C14) 



Therefore we find that the wave function in the presence of the single impurity at a; = x^ remains unchanged from 
the original one 



4>n{x) 













+ 


_^^-i(p ^ikx 



(C15) 



The first term represents the incoming wave at the electron branch. The second term expresses the outgoing wave 
in the hole branch. The Andreev reflection is perfect and the normal reflection is absent even in the presence of the 
single impurity in the normal metal. With using the Blonder-Tinkham-Klapwijk formula, 



G 



ATS 



[1 



„ee 1 2 



„he |21 



(C16) 



the zero-bias conductance of the NS junction remains unchanged from Gns = jh independent of the impurity 
scattering. 

On the way to the conclusion, we derive a relation 



eh he _ i 
' nn nn ^ 



(C17) 



at _E = 0. This plays an important role in the resonant transmission of a quasiparticle in a normal metal. For 
comparison, the Andreev reflection coefficients of the s-wave transparent NS junction becomes 



at = 0. However, they gives a relation r; 



„he 



ih „/ie 



„eh 



-le 



lip 



(C18) 



= — 1. In this case, a quasiparticle is scattered by the impurity and 



the conductance decreases from G^s — Gq. The relation in Eq. (|C17p is equivalent to the necessary condition for 
the formation of Andreev (Majorana) bound states at = 0. 
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Appendix D: Analysis of quasiclassical Usadel Equation 

In this section, we consider a diffusive normal metal is attached to p^-wave superconductor as shown in Fig. 1(d). 
At first, we define the quasicalssical Green functions in terms of Gor'kov Green functions. In the mixed representation, 
Gor'kov Green functions become 

Gix,t;x' ,t') ^Gixcx - x' ,tc,t ^ t') = J gG(xc, fc, t^, e)e*'=(^-^')-^^(*-*'), (Dl) 

Fix,t;x',t')=Fix,,x-x',U,t-t')^ J^J ft, i,, e)e^'=(---')— (*"*'), (D2) 

x + x' t + t' 

= ^ , = — (D3) 

When we consider static state, the Green functions are independent of tc- With replacing Xc by x, the quasiclassical 
Green functions are defined as 

g(x,fc,e)=- f d^kG{x,k,e)~- [ d^k^, (D4) 
fix,k,e)^- f d^kF{x,k,e). (D5) 



They obey so called Eilengerger equation. In what follows, we fix the phase of the superconductor ip at 0. When the 
normal metal is in the dirty limit, g{x, fc, e) and f{x, k, e) are isotropic in momentum space. Since they satisfy the 
normalization condition g^{x,e) + p{x,e) — 1, it is possible to apply a parameterization: g{x,e) — cos[6'(2:, e)] and 
/(x, e) ~ sin[0(a;, e)]. The function 6{x^ e) obeys the Usadel equation in the diffusive normal metal 

D ' +2iesin[e{x,e)]=0, (D6) 

where D is the diffusion constant in the dirty normal metal. 

In what follows, we consider NS junction shown in Fig. 1(d), where a dirty normal metal is introduced between 
a clean normal lead wire (x < —L) and a pj;-wave superconductor (x > 0). The boundary condition for 0{x,e) are 
given by2^ 

0{x = ~L,e) ^0, (D7) 



L d9{x,e) 



Rn dx 



fs cos 6*0 - gs sin 6*0 



x=0 



Rb 2 - |i„|2 + |i„|2(/5sin0o +5scos^o)' 



with 



a ai n ^ 9+ +9- . ■ f+9- - f-9+ /t^^n 

=6* X = 0, e), gs = — ^ „ „ , fs = ^— . , , , (D9) 

^ + 9+9- + f+f- ^ + 9+9- + f+f- 

Rb = [GqK\''Y\ ^q^\- (DIO) 

The parameter Rf^ and Rb are the normal resistance of the dirty normal metal and that due to the potential barrier 
at the NS interface, respectively. The information of the pairing symmetry of superconductor is embedded in the 
surface Green function gs and fs- The total resistance of the junction R is represented byS^ 

R=Rb + Rn, (Dll) 
1 Gq 

"2 1(2 - |t„|2) + \tr,\Hcoseogs + sinOofsW ' ^^^^^ 
Co =|i«P(l + I cos flop + I sinflon(l + \9s\^ + \fs\^) 

+4(2 - |<„n [Rc(5s)Re(cosflo) +Re(/s)Re(sinflo)] +4|i„pIm(cosflosin*flo)Im(55/5), (D13) 

^Rn f° 2dx 

^ L J_Ll + \cos0{x,e)\^ + \sm0{x,e)\^' ^ ' 

The resistance Rb and R^ are not equal to their normal one's Rb and Rj^. They are modified by the proximity 
effect. 
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In the case of the px-wave superconductor, following relations hold 

ff+ = 5- = ^ , /+ = -/- = . (D15) 

The Pa;-wave symmetry of superconductor is represented by the relation /+ = —/-• At the surface of Px-w&ve 
superconductor, purely odd-frequency pairing state exist due to the formation of Andreev bound state as discussed 
in Sec. 2.3. 
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FIG. 6: Conductance are plotted as a function of L in NS junctions of nano wire. 



The Usadel equation can be solved analytically at e = 0. Under the boundary condition {L / RN){d6 / dx)\x=a = iGq, 
we obtain 



e{x,e = 0) ^iRpfGq 



x + L 



]{x) —cos6{x,0) = cosh ( RnGq 



L 



f(x) =sin6'(a;, 0) — isinh ( RnGq 



L 

x + L 



(D16) 
(D17) 

(D18) 



at e = 0. The pairing function f{x) represents the spin-triplet s-wave odd-frequency pair. Indeed, f{x) is purely 
imaginary number at e = 0. Finally we obtain the zero-bias resistance 



R —Rb + Rn — Gq , 

Rb =7^ [1 + »/(0)/.g(0)] = 7^[1 - tanh(GQi?jv)], 
Gq Oq 

Rn =7^ [-if[Q)/gm - tanh(GQi?^). 
Gq Gq 



(D19) 
(D20) 

(D21) 



The total resistance at the zero-bias voltage is independent oi Rn and Rb, and remains unchanged from R = Gq . 

It is worth to consider the physical meaning of the resulting resistances Rb and Rn- Rb is the resistance at the 
interface which decreases from Gg^ with the increase oi Rn- In other words, the interface resistance decreases with 

the increase of the amplitude of odd-frequency pair /(O). Rn is the resistance of the dirty normal metal. In the limit 
of weak proximity effect RnGq <C 1, the amplitude of odd- frequency pairs becomes small. In such limit, we find 
Rn = Rn- On the other hand for RnGq ^ 1, —if{0)/g{0) goes to unity and Rn approaches to Gq^. Thus the 

odd-frequency pairs play a crucial role in the relation of i? = Gq^. 

The diffusive transport is assumed in the Usadel equation. It is possible to check the validity of Gns — Gq when 
the normal segment is in the localization regime. In Fig. [51 we plot the zero-bias conductance in NS junctions of nano 
wire as a function of the length of the disordered segment L. When the nano wire is non topological at A = Vex = 0, 
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the zero-bias conductance decreases exponentially with L due to the localization. On the other hand, the conductance 
for topological nano wire junctions at A = 0.5i and Vex — l-5i remains unchanged from Gq even in the localization 
regime. 
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